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Abstract. 

Considering a multi-pathway structure in a light-harvesting complex of 
photosynthesis, we investigate the role of energy- level mismatches between antenna 
molecules in transferring the absorbed energy to a reaction center. We find a condition 
in which the antenna molecules faithfully play their roles: Their effective absorption 
ratios are larger than those of the receiver molecule directly coupled to the reaction 
center. In the absence of energy- level mismatches and dephasing noise, there arises 
quantum destructive interference between multiple paths that restricts the energy 
transfer. On the other hand, the destructive interference diminishes as asymmetrically 
biasing the energy- level mismatches and/or introducing quantum noise of dephasing for 
the antenna molecules, so that the transfer efficiency is greatly enhanced to near unity. 
Remarkably, the near-unity efficiency can be achieved at a wide range of asymmetric 
energy- level mismatches. Temporal characteristics are also optimized at the energy- 
level mismatches where the transfer efficiency is near unity. We discuss these effects, 
in particular, for the Fenna-Matthews- Olson complex. 
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1. Introduction 

Photosynthetic complexes are sophisticated light-harvesting machinery consisting of 
antenna molecules. The energy absorbed by the antenna molecules, so-called exciton, is 
transferred to many intermediate molecules and eventually arrives at a reaction center 
(RC) where the process of biochemical energy conversion is initiated. Recent work 
has reported that quantum theory governs the exciton transfer in some light-harvesting 
complexes that harness the absorbed energy with almost 100% efficiency [I]-[3]. It 
has also been suggested that the interplay of quantum walk and noisy environment 
provides the high light-harvesting efficiency in these complexes [4]-[T5]. However, 
the underlying mechanism of noise-assisted enhancement has remained elusive and 
the role of structural characteristics of antenna complexes is still not clear. First, 
photosynthetic complexes possess energy-level mismatches between antenna molecules. 
As energy-level mismatches are likely to cause Anderson localization [16] and to 
inhibit the transfer of excitation it is desired to understand why photosynthetic 
complexes evolve by maintaining such energy-level mismatches instead of eliminating or 
reducing if impossible. Second, the noise-assisted enhancement is not possible in certain 
situations, such as uniform linear chains with no energy-level mismatches [5j. These 
observations suggest that energy-level mismatches are related to the condition for the 
enhancement [6], but their fundamental role is still not unraveled. Moreover, even in 
the cases where the noise-assisted enhancement is possible, the underlying principle of 
enhancement remains ambiguous due to the lack of basic studies based on quantum 
interference [7j. 

In this paper, we investigate the role of energy-level mismatches in a multi-pathway 
complex where multiple sub-complexes are independently connected to the RC via a 
receiver molecule (see figure [T]). For a single-pathway complex, we show that any energy- 
level mismatches suppress the efficiency of exciton transfer as expected by Anderson 
localization. In the absence of energy-level mismatches, we show that quantum noise 
of dephasing, which decreases the quantum coherence of excitation, will never help 
the transfer efficiency. For a bi-pathway complex, on the other hand, we demonstrate 
that the absence of energy-level mismatches leads to undesired destructive interference 
at the receiver molecule. The destructive interference is caused by the two types of 
probability amplitudes, one of coming from one pathway and the other of going to and 
returning from the other pathway. It blocks exciton transfer to the RC so that the 
transfer efficiency is less than 50% and relaxation of the sub-complexes to the ground 
state becomes very slow. We show that energy-level mismatches play a dominant role in 
suppressing the undesired interference. Then, the presence of energy-level mismatches 
enhances the transfer efficiency even though the resultant localization effect suppresses 
the energy transfer in each pathway of the complex. Due to the competition between 
the Anderson localization and the destructive interference at the receiver molecule, a 
moderate amount of energy-level mismatches will improve the light-harvesting efficiency 
as minimizing both the undesired effects in the energy transfer. This picture is consistent 



Role of Energy-Level Mismatches in a Multi-Pathway Complex of Photosynthesis 



3 



with our results. In addition, dephasing noise is found to relaxe the localization and to 
suppress the undesired interference at the receiver molecule so that the cooperation 
of energy-level mismatches and dephasing noise significantly improves the transfer 
efficiency after all. 



2. Light-Harvesting complex 

Light-harvesting complex is modeled as a system consisting of n two-level molecules, 
whose dynamics is governed by a master equation in the form of 

j t P = ~[H, p] + L A (p) + L D (p) + Ldp(p), (1) 

where p is the density matrix of the molecules and H is the Hamiltonian of the system, 
given by 

n n 

H = J2 m j( T+aJ + H"fc«^fc + °t<rj)- (2) 

j=l j<k 

Here, and a J are the raising and lowering operator for molecule j, Mlj is the excited 
energy of the molecule and Jjk is the electronic coupling constant between molecules j 
and k. The first and second non-unitary terms La(p) and L D (p) describe the processes 
of absorbing and emitting thermal light and phonons with coupling constants rjj and 
Tj at molecule j. The second term L D (p) also contains an irreversible decay from the 
receiver (denoted by j = 1) to the RC with a coupling constant T R c. The Lindblad 
operators are given by 

n 1 

L D (p) ^Mvjpvf - \i<rI<Tj,p}), (4) 
j=i 

where aj = Nirjj + N p Tj and /3j = (Ni + Vjrjj + (N p + l)Tj + T R c^ij are the exciton- 

creation and exciton-decay constants of molecule j. Here, JVj (N p ) is the mean photon 

(phonon) number of the thermal light (phonons) and Sij is the Kronecker delta with 5\j 

indicating that the RC is coupled with molecule 1. The last non-unitary term L DP (p) 

describes the dephasing process, given by 

n 1 

^dp(p) = ^27j^t a 7P a t a 7 ~ 2^ a J +a J 7 '^)' ( 5 ) 

where jj is the dephasing constant of molecule j. This process is also caused by the 
interaction with phonons with conserving the system energy. By these non-unitary 
processes, the system state decoheres. To investigate continuous exciton transfer, we 
consider the steady state of the master equation ([T]), which is invariant as time evolves. 
To clarify the principal mechanism, we begin with a single-pathway complex consisting 
of two antenna molecules and then consider a bi-pathway complex of three antenna 
molecules (see figure [j]) . 
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Figure 1. Schematic representations of (a) the single-pathway complex and (b) bi- 
pathway complex. The molecule coupled with the RC is called a receiver and the 
other ones donors. The solid lines between molecules represent electronic coupling 
that enables an exciton to hop between them. 



3. Single- Pathway Complex 

3.1. Steady Exciton Transfer 

We shall express the steady state condition -^p = in terms of the probabilities by 
eliminating the coherence terms. For the single-pathway complex in figure [jja), the 
steady state is given in the form of 

p = P |0> <0| + Pi |1) <1| + P 2 12) <2| + Ci 2 |1) <2| + C? 2 |2) (1| + P 12 1 12) (12| , (6) 

in localized exciton basis {|0) , \j) = |0) , 1 12) = afa^ |0)}, where |0) is the ground 
state of the system. Through the master equation ([I]), the steady state condition jj^p = 
results in the set of equations 
dP 



dt 
dP 1 

~dT " 
dP 2 

dt 

dCu 

dt 
dC* l2 

dt 
dP 12 



-(a 1 + a 2 )P + p 1 P 1 + p 2 P 2 = 0, (7) 

«iPo - («2 + Pi)Pi + iJu{C l2 - C* 2 ) + /3 2 Pi2 = 0, (8) 

« 2 Po - («i + /? 2 )P 2 - iJi2(Ci2 - C{ 2 ) + AP12 = 0, (9) 

= iMPi - P 2 ) - (P12 + *Aft 12 )C 12 = 0, (10) 

= -i Ji 2 (Pi - P 2 ) - (P12 - iA0 12 )C* 2 = 0, (11) 

= a 2 P x + a x P 2 - (A + /3 2 )P 12 = 0, (12) 



dt 

where Pi 2 = \{ct\ + fix + 71 + a 2 + f3 2 + 72) and Af2i 2 = Oi — f2 2 . The equations for the 



coherence terms, equations (10) and (11), lead to 
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By substituting equations (13) and (14) into equations (JsJ) and Q, the steady state 
condition 4-p = can be expressed in terms of probabilities 



dtf 

dPo 

dt 
dP 1 



dt 
dt 

dPu 
dt 



-(a 1 + a 2 )P + APi + &P2 = 0, (15) 

aiP - (« 2 + f3 1 )P 1 + 62(^2 - Pi) + /? 2 Pi2 = 0, (16) 

« 2 Po - («i + &)P 2 + 6 2 (Pi - P 2 ) + AP12 = 0, (17) 

= a 2 Pi + «iP 2 - (A + /3 2 )Pi2 = 0, (18) 

where the hopping constant £12 is defined by £12 = Ao2 ■ Focused on ^Pj, c^-Po 

(fSjPj) is the exciton-increasing (decreasing) rate of Pj from (to) the ground state, and 
PkPi2 ( a kPj) is the converting rate of Pj from (to) the two-exciton state. In addition, 
^uPk is the exciton hopping rate from molecule k to molecule j, while £12 Pj is that 
from molecule j to molecule k. It is notable that the hopping constant £ 12 contains 
significant information, such as electronic coupling constant J12, energy-level mismatch 
between the donor and receiver A^i 2 and dephasing constants 7^. The solution to 



equations ( 15)— ( 18) is given, with the normalization Tr(p) = Pq + Pi + P2 + P12 = 1, as 



PMm + a 2 + fa + fa) + {fa + fafi 



Pv = 7 o^Ttu ^ ' ^7 ' ^ 2 5 ZP^p ( 19 ) 



(ai + a 2 + /?i + /3 2 )[(c*i + /?i)(a 2 + AO + («i + «2 + A + fa)^} 
p aifafai + Qfr + fti + ft 2 ) + (gi + a 2 ){fa + /3 2 )Ci2 

1 " (ax + a 2 + /?i + /5 2 )[(«i + fa){a 2 + A) + («i + OL2 + P1 + fa)£u] ' 
p _ Qfrftifoi + ^2 + fti + ft 2 ) + (gi + a 2 ){fa + /3 2 )Ci2 

2 ~ (ai + a 2 + fa + /3 2 )[(«i + Pi){ol 2 + fa) + (ai + a 2 + A + ^2)62] ' 
p _ ^1^2(^1 + a 2 + /3i + fa) + (ai + « 2 ) 2 Ci2 

12 " (c*i + a 2 + /3i + /3 2 )[(c*i + /3i)(«2 + A) + (c*i + a 2 + /3i + /3 2 )fi 2 ] 

Necessary and sufficient condition for the faithful donor 



We shall find the condition in which the antenna molecules faithfully play their roles. 
An exciton is created at molecule j, if unoccupied, by the rate of E$ = aj(l — Wj) 
where Wj = Pj + P\ 2 is the probability of finding an exciton at molecule j, given by 

w = cx 1 (a 2 + fa) + {a 1 + a 2 )£ 12 . . 

1 {a x + fa){a 2 + fa) + (c*i + a 2 + fa + /3 2 )fi2 ' 

w = + (^1 + ^2)62 , 24 x 

2 (ax + fa){a 2 + A) + {on + a 2 + fa+ fa)£ 12 

On the other hand, molecule j loses its exciton, if occupied, by the rate of = PjWj. 
Here, R$ includes the transfer to the RC that happens by the rate of i?Rc = TrcWi. 
The net rate of exciton hopping from the donor to the receiver R^ = £12(^2 — Pi) is 
given by 

^(12) (® 2 fa ~ OLlfa)£l2 
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so that we have two rate equations: i?« + = R$> and = R$ + R% 2 \ 

Now we consider the effect of attaching the donor to the receiver (£ 12 > or 
equivalently J 12 7^ 0) compared to the decoupled system (J 12 = 0). The attachment 
increases the net absorption rate of light (and also phonons) Ra = + Ra\ as seen 
from the positivity of the derivative of the net absorption rate Ra with respect to the 
hopping constant £12 



dR, 



12 



a 2 fa — ai/3 2 



2 



Xa 1 + fa)(a 2 + fa) + (a 1 + a 2 + fa + fa)£ 12 \ >0 ' ^ 
One might conjecture that the attachment also increases the exciton transfer rate to 
the RC, but this is not necessarily the case. In fact, the transfer rate is enhanced if and 
only if 

dR R c r RC (^2 + fa)(a 2 fa - aifa) , 2? x 

d£i2 ~~ + fa)(a 2 + fa) + («i + ^ 2 + A + /3 2 )£ 12 ] 2 ' 1 j 

or equivalently 

A 2 > A u (28) 

where Aj = &j//3j is the effective absorption ratio of molecule j. This condition holds 
if T R c is large enough as A\ is proportional to T^L This is also the condition that 
the net exciton hopping is directed from the donor to the receiver {R% 2) > 0). Thus, 



attaching a donor in the condition (28) faithfully directs the energy flow from the antenna 



molecules eventually to the RC as well as increasing the net absorption rate of light. If 



the condition (28) is not satisfied, the net exciton hopping is directed from the receiver 
to the donor (w H < 0) so that removing a donor from the antenna complex ( J 12 = 0) 
increases the exciton transfer rate to the RC. These results are clearly reflected in the 
transfer efficiency e = Rrc/Ra & [0, 1], given as 



total 



, fa 0^2 fa ~ Oil fa 



(29) 



Aotal £12(^1 + Ot 2 ) + ^l(^2 + fa) 

where /3 to tai = Y^=i Pj ^ s the total exciton-decay constant of the antenna complex. 

We would note that the attachment decreases the receiver's absorption rate R^ while 

(2) 

increasing the donor's R A \ thus reducing the role of receiver in the absorption and 
changing it to the transmission to the RC. Under physiological conditions, the intensity 
of sunlight is weak (aj <C fa) and the single-exciton manifold is of primary importance 
for modeling photosynthetic complexes such as the Fenna-Matthews-Olson (FMO) 
complex [T71 HE]. In this case, the steady state is well-approximated within the single- 
exciton manifold and e is reduced to the single-exciton transfer efficiency ji]-[T5]. 



POO 

Trc / dt(l\p(t)\l), 
Jo 



(30) 



where the initial state is given by p(0) = Z)j=i( a jMotai) \ j) 01 with a total = Y^=i a j 
and its dynamics is governed by the master equation in the form of 

j t P = - l -[H,p]+L D (p) + L DP (p). (31) 
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Figure 2. (color online) Transfer efficiency e of figure [TJ a) under low energy absorption 
(<*i = 0, a 2 < /? 2 ). We choose r RC /Ji 2 = 10" 1 , T 3 /J 12 = 1(T 3 (fa = Tj + r RC ^), 
I3/J12 = 1/J12 an d A^i 2 /Ji2 = k for k = 0, 1, • • • , 5. The region of large 7/J12 is 
only for eye-guiding. 



3.3. Subsidiary Role of Dephasing Process 

We will consider the noise effect on the transfer efficiency e with respect to the degree 
of the energy-level mismatch, as shown in figure [2j Assume that the condition of 



enhancement for the transfer efficiency, Eq. (28) holds. Adding dephasing noise increases 



e when the energy-level mismatch |A^i 2 | is larger than the constant D 12 of no dephasing 
7j = 0. However, e is still maximized when there are no energy-level mismatch and no 
dephasing (AQ12 = 7j = 0). The noise, if present, never makes any enhancement. 

These results remain unaltered for a single-pathway complex consisting of more 
than two molecules. Plenio et al. [5] investigated the transfer efficiency for a single- 
pathway complex of a uniform linear chain where all molecules have the same energy 
level, and their electronic couplings and exciton-decay rates are uniform. Assuming 
that the initial excitation is located at an edge of the chain and the receiver molecule 
is located at the opposite one, they numerically observed that the dephasing does not 
improve the transfer efficiency. The observation implies that there is no Anderson 
localization effect and the dephasing noise never enhances the transfer efficiency if any 
energy-level mismatches are absent. Large energy-level mismatches cause the Anderson 
localization and the dephasing noise releases the inhibition due to the localization, 
leading to the enhancement of the transfer efficiency. However, the optimal structure 
of a single-pathway complex is that with no energy-level mismatches. In this sense, the 
noise-assisted enhancement of e is a subsidiary effect once the energy-level mismatch 
settles in the system. In the next section, we show that the negative role of energy- 
level mismatches and the subsidiary role of dephasing noise dramatically change in a 
multi-pathway complex. 

4. Bi-Pathway Complex 

4-1- Enhancement of Transfer Efficiency by Energy- Level Mismatches 

Consider a bi-pathway complex as in figure [jjb), where two identical donors (o^ = a 2 = 
<^3, Pd — 02 = 03: la = 72 = 73) are independently coupled to a receiver molecule 
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(J23 = 0). We assume low energy absorption (ctj <C taking the single-exciton 
approximation. We also take a\ = 0, reminding that attaching donors transforms the 
receiver's role to a transmission channel in the single-pathway complex [|} In this case, 
the transfer efficiency e is reduced to the single-exciton transfer efficiency when the 
initial state is given by |(|2) (2| + |3) (3|): 



RC 



total 



1 + 



PlPd §(£l2 + £l 3 ) +2(63 + I A) 



-1 



(32) 



Aotal ^12^13 + (£l2 + £13) (£23 + \Pd)_ 

where /?totai = Pi + P2 + 03- Here, the hopping constants ^ between molecules j and k 
are given by 

t ^12(^12 _ ^13 + ^13523) /qq\ 



J 12^ 21 + ^13^13 + S21S13S23 

^13(^13 ^12 + ^12^32) 
^13^31 + ^12^12 + S31S12S32 



t _ ^13W13~ ^12 ^^12^32; 



^23 = 72 C,72Q,CCC + C ' C -' ( 35 ) 
^12^21 T ^13^13 1 ^21^13^23 

where c.c. stands for the complex conjugate, Sjk = Djk+iAQjk, Djk = \{Pj+lj+Pk+lk) 
and AQjk = Qj — Q^. Note that e is a monotonically increasing function of the hopping 
de 

constants ^— - > 0, Vj, k. 

The transfer efficiency e is presented as a function of energy-level mismatches in 
figure [3| Note that in each panel, near-unit transfer efficiency is located at a wide range 
of energy-level mismatches and the maximal transfer efficiency is located on asymmetric 
energy-level mismatches AQ12 7^ A^i 3 . In the first column of panels where no dephasing 
is assumed, the effect of energy-level mismatches are clearly shown: Symmetric energy- 
level mismatches A^i 2 = A^i 3 result in the suppression of e < 1/2, whereas asymmetric 
ones lead to the near-unit transfer efficiency. As in the second column of panels, adding 
dephasing 71 at the receiver can slightly enhance e, but the previous tendency by the 
energy-level mismatches is rather unaltered: For the symmetric energy-level mismatches 
AQ = AQ12 = Af]i3, in the absence of dephasing noise at the donors (7^ = 0), the 
transfer efficiency is reduced to 



RC 



Pi+Pd 



< \ (36) 



where Did = |C#i + 71 + Az). However, it changes dramatically when the donors 
are under dephasing noise. Increasing dephasing 7^ at the donors results in the great 
improvement so that e goes over 1/2 near to the unity even for the symmetric energy-level 
mismatches. These results remain valid for asymmetric coupling constants «7 12 7^ J13, 
as seen by comparing three rows in figure [3j 

X In both cases of figure[lja) and (b), under the same condition (28), the net exciton transfer is directed 
from the donor (s) to the receiver. 
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Figure 3. (color online) Transfer efficiency e of the bi-pathway complex in figure [TJb), 
parameterized by energy-level mismatches. In each panel, horizontal and vertical axes 
represent energy-level mismatches A^i2/Ji3 and Afii3/Ji3, respectively. From first 
to third rows, we take J12/J13 = 1,2,4, respectively. In both cases, the parameters 
are chosen as = 10 _1 , Tj/Jis = 10 -3 (/3j = Tj + Trc^ij) an d the dephasing 

constants (71/ J 13 , j d / J 13 ) = (0,0), (10 _1 ,0), (0, 10" 1 ) from the left to the right. 



4.2. Destructive Interference at the Receiver Molecule 



The underlying mechanism is the quantum destructive interference at the receiver 
molecule and its suppression by the energy-level mismatches and/or the dephasing noise. 
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In the presence of energy-decay process, once excited by light absorption, a localized 
state \j) evolves into a statistical mixture, (1 — Pj(t)) |0) (0| + Pj(t)pj(t), of the ground 
state and a delocalized state Pj(t). Here, the delocalized state is a superposition of 
localized single-exciton states. In the absence of dephasing noise, the delocalized state 
Pj(t) remains in a pure state p jit) p jit) = \^jit)) {^j{t)\ f° r which it is convenient to 
employ a stochastic Schrodinger equation [19l [20] , 

ihj t \^j(t)) = K\^(t)) with K = H-ihJ2\^t^. (37) 

k=l 

Here, \i/jj(t)) is unnormalized and Pjit) = {^jit)\^jit)) is the probability that an 
exciton is still present in the molecules. For the symmetric energy-level mismatches 
AQ = A^i2 = Af2i3, the delocalized state is written as 

i*w> = ^«p(§ t )i«.>+E^°p<§')i«v>. < 38 > 

p=2 F 

\Ut)) = exptf*) M + £ ^ exp(| t ) W , (39) 

where |^(0)) = \j) for j = 2, 3. Here, \v p ) are the right eigenvectors of given as 

|vi> = Ji 3 |2> - Ji 2 |3) , (40) 
A P 

K) = + J 12 |2) + J 13 |3), p^l, (41) 

where AEi p — E\ — E p , M p = J\ 2 + Jf 3 + (AEi p /h) 2 and £" p are the complex eigenvalues 
corresponding to \v p ), given by 

E 1 =ih(-i(Q 1 -AQ)-^ d ), (42) 
E p = ih(X-(-l)>Y), p^l, (43) 



where X = - \AQ) - \{fa + fa), Y = \^ - /3 d + 2*Af2) 2 - 1Q(J 2 12 + J? 3 ). 

Note that when normalized, has no probability amplitude of staying in the receiver 
due to the perfect destructive interference between the amplitudes of multiple paths to 
and from the donors. Thus, does not contribute to the energy transfer to the RC, 
so that it states the eventual loss of excitons at the donors. The destructive interference 
consequently results in the transfer efficiency less than 1/2 (see figure [i]). For the 
asymmetric energy-level mismatches, on the other hand, the destructive interference 
diminishes as every eigenvector of K has amplitude at the receiver. The transfer 
efficiency increases as biasing the asymmetric degree of energy-level mismatches and 
turns to decrease for large mismatches (see figure [3]) . The destructive interference can 
also be avoided by the dephasing at the donors that collapses probabilistically into 
one of the localized states at the donors, which have a chance to transfer an exciton to 
the RC. 
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Figure 4. Population relaxations of molecules 1, 2 and 3 with respect to time in 
the bi-pathway complex when there are no energy-level mismatches and no dephasing 
noise, but the decay constants Trc/J = 10 _1 and Tj/J = 10 -3 (/3j = Tj + Trc^ij). 
An exciton created at molecule 2 hops over the receiver and donor molecules, and 
eventually relaxes in two different ways: Either the exciton is transferred through 
the molecule 1 (receiver) to the RC, or it relaxes at the donors to the environment. 

(a) The population at the receiver molecule becomes negligible for Jt > 60, whereas 

(b) /(c) the populations at the donors last even beyond 60 units of Jt. These results 
imply the evidence of the destructive interference at the receiver molecule. Due to 
the destructive interference, since about 60 units of Jt, the exciton is trapped at the 
donors and consequently relaxes to the environment. 



4-3. Quantitative Description of Noise-Assisted Enhancement 

Quantitative description of the noise-assisted enhancement can also be provided by using 
a stochastic Schrodinger equation. Assume a localized state \j) is created at time t — 0. 
In the presence of the exciton-decay and dephasing processes, the localized state \j) 
evolves into a statistical mixture of the ground state and a delocalized state, 

p=(l-p i (t))|0>(0|+p i (t)p i (t). (44) 

Here, (1 —Pj(t)) and Pj(t) are the probabilities for finding the system in the ground state 
and the delocalized state, respectively. The delocalized state Pj(t) is a statistical mixture 
of single-exciton states. As time evolves, the delocalized state Pj(t) is probabilistically 
collapsed into the ground state |0) by the exciton-decay process at molecule k with 
the rate of fa (k \ p(t) \k). The delocalized state Pj(t) is also probabilistically collapsed 
into the localized state \k) by the dephasing process at molecule k with the rate of 



lk(k\p(t) \k). These stochastic processes are equivalent to the master equation (31) 
that can be rearranged as 

d i n n 

J t P=-^( K P- P R] ) + ^2^kP (7 k + ^Zlk(rt (J kP (J t (J k^ (45) 

k=l k=l 

where K = H — ifoY^k=i ^Wk + 7k)^k a k • To describe the dynamics of the state p based 
on the stochastic processes, we express the delocalized state pj(t) as a statistical mixture 
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of pure single-exciton states 

n pt 

Pj (t) Pj (t) = WMiMtM + Y, / ^il^(*;*i)X^(*;*i)I^S P) (*i;0) 



k=l 



k,i=i Jt ± Jo 

Here, \ipk(t;t z )) is an unnormalized pure state such that \^k(t]t z )) = for time t < t z , 
\^k(t;t z )) = \k) at t — t z and the dynamics of \ipk{t]t z )) for t > t z is determined by a 
stochastic Schrodinger equation [19j [20] in the form of 

ih^\M*Q) = K\MWz))- ( 47 ) 

The formal solution is given as \^k(t]t z )) = exp[— |K(£ — t z )] \k) = \ipk(t — t z ]0)) for 
t > t z . R^ F \t z+ i]t z ) = 7/ \(l\i/; k (t z+1 ]t z ))\ 2 is the rate that \ijj k (t]t z )) is collapsed into 
the localized state \l) by the dephasing process at t = t z +\. The probability to be in the 
ground state |0) at time t is given by 



n pi n pt pt 



(DP) 



(*i;0) 



n 

+ 53 / efts / dt 2 / <MS^s;*2^^ (48) 

k,l,m=l Jt * Jt! J0 

where R^\t z+ i]t z ) = f$\ \ (l\ipk(t z +i]t z ))\ 2 is the rate that |^fc(t; £ z )) is collapsed into the 
ground state by the exciton-decay process at t — t z +\. 

We shall express the single-exciton transfer efficiency e in equation (32) as a series 
sum of the transfer probabilities to the RC on the basis of the stochastic processes. There 
happen two types of stochastic processes in our model, one is the exciton-decay process 
and the other is the dephasing process. In order to analyze what processes and how many 
times the processes are involved in the transfer to the RC, it is convenient to consider 
the temporally- accumulated stochastic probabilities. A (temporally- accumulated) decay 
probability is defined by 

POO POO 

pi?= mmi>iftt z ))\ 2 = / dmm^it-m 2 v*, (49) 

Jt z Jo 
where we used \^ k (t]t z )) = \iftk(t ~ t z ]fy) for t > t z . This is the transition probability 
from the state \ j) to the ground state |0) conditioned by a decay process at molecule k. 
Similarly, define a (temporally-accumulated) dephasing probability as 

POO POO 

P tT )= dt lk \{k\^{t-t z ))\ 2 = dt lk \{k\^(t-Q))\ 2 yz, (50) 
Jt z Jo 

which is the transition probability from the state \j) to the state \k) conditioned by 

a dephasing process at molecule k. Here, includes both probabilities of the loss 

to the environment and of the transfer to the RC through the receiver molecule. The 

transfer probability to the RC is segreg ated by {T KC /Pi)P[f \ As there are no other 
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stochastic processes assumed, the probabilities P^p and p£? P ^ satisfy the normalization 



condition, 



kj 



where P 1 



(S) 



k=i 



kj 



P 



kj 



P, 



(DP) 



kj 



(51) 



f™dtD kk \(k\^j(t]0))\ 2 and D kk = f3 k + 7*.. Now, the 



single-exciton transfer efficiency 6j of the initial state \j) is expanded into the sum of 
the stochastic chains of the decay and dephasing processes: 



(r RC /ft) 



E A<f 4 DP) + E < , 4 DP, ^ DP| + - ( 52 > 



k = l 



k,l=l 



(r RC //3i) VFd(Pdp) 2 



,2=0 



(53) 



Here, a single chain (r R c//3i)(PD(PDp) 2: )ij is the transfer probability to the RC of the 
initial state \j) through z times dephasing processes followed by the decay process. The 
matrices P& and Pdp of the decay and dephasing probabilities are given, respectively, 
as 



Pn = 



-Pdp — 



Pi 



7i 







Pi 




72 








P 3 

0" 



73 



-D u - 62 
62 
63 

■Du - 62 ■ 

62 
63 



13 



;i3 



62 
-D22 — C12 

<^23 

62 

-D 2 2 — £l2 

^23 



^23 



63 

^23 

— ^33 £l3 £23 



-1 -1 



^23 



63 

£23 

— ^33 £l3 ^23 



-1 -1 



(55) 



where ^ are the hopping constants defined in equations (33)-(35). The single-exciton 



transfer efficiency e in equation (32) is expressed in the form of 

3 

e = E' 



1 



tj x -. 



(56) 



where 1/2 is the probability that a single exciton is initially created at molecule j. 

The stochastic chains can be applied to explain the noise- assisted enhancement. 
Particularly, consider the case of symmetric energy-level mismatches AQ = AQ12 = 
The sum of the probabilities that the initial exciton decoheres to the donors 
(molecule k — 2, 3) once created at the donors (molecule j = 2, 3) in the probabilities 
1/2 is given as 



p 



(S) _ 



dd 



En, 

j,k=2 



(S) 



X 2 = 1 



2 + D. 



dd 



+ 



Du + Atf 



Du Jl 



12 



J 2 
J 13 



D 



id 



1 

> 2' 



(57) 



where D ld = \{pi + 71 + d + Id) and D dd = f3 d + ld . Noting P fc ( f = (f3 d /D dd )P^ and 

( S) 

for k = 2, 3, the decay and dephasing probabilities at the donors 



p(DP) 
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Table 1. Transfer efficiency and temporal characteristics of the bi-pathway complex 
used in figure [5] and figure [6] for no energy- level mismatches AQ12/J = AQ13/J = 
and anti-symmetric energy-level mismatches AQ12/J = — AQ13/J = ±1. 



Afl 12 /J 








(trcJ)" 1 


(trJ)- 1 


e 














0.051 


0.002 


0.490 








10- 1 





0.051 


0.002 


0.490 











10- 1 


0.026 


0.026 


0.962 


±1 


Tl 








0.034 


0.034 


0.971 


±1 


Tl 


10- 1 





0.034 


0.034 


0.971 


±1 


Tl 





10- 1 


0.034 


0.034 


0.971 



are respectively given as 



3 



p(D) = p (D) 1 _ _^_ p (S) p (DP) p (DP) 1 _ J^ p (S) /roN 

r dd — r kj o n dd ' dd ^' 9 — n dd ? w ; 

J,k=2 j,/c=2 

while those at the receiver molecule are respectively given as 



3 3 



Provided there is no dephasing noise at the donors (7^ = 0), the decay probability at 
the donors pffi = P^ and it is larger than 1/2. Reciprocally, the decay probability at 
the receiver is less than 1/2 and so is the transfer probability to the RC. This explains 



the transfer efficiency e < 1/2 in equation (36). On the other hand, if introducing and 



increasing the dephasing noise at the donors (7^ > 0), the decay probability at the 
donors decreases to zero, pffi = (fid/ D^P^ —> 0, whereas the dephasing probability 
at the donors P^^ increases. Then, as seen in equation (52) or (53), the stochastic 



processes are summed up to result in the high transfer efficiency. 

4-4- Temporal characteristics of Relaxation to the Ground State and of Transfer to the 
Reaction Center 

We shall investigate temporal characteristics of the single-exciton transfer. In order to 
consider relaxation process that the initially excited state decays to the ground state, 
we define relaxation time tr as 

poo n 

TR= / dtTPj(j\p(t)\j)t. (60) 

Jo j=1 

Here, Y^j=\Pj 01 \ j) ^ s the probability density function of time t that the system 
state p(t) decays to the ground state at i, satisfying the normalization condition, 
J °° dt YTj=i Pj 01 \ j) = 1- ^ the probability of finding an exciton in the system 
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(c) Dephasing at the donors 



z zz 

0.026 0.034 0.051 

Figure 5. (color online) Inverse transfer time to the RC (trcJ) -1 of the bi- 
pathway complex in figure [TJb) with symmetric couplings constants J = J12 = «/i3, 
parameterized by energy-level mismatches. In each panel, horizontal and vertical axes 
represent energy- level mismatches AQ12/J and AO13/J, respectively. The parameters 
are chosen as T^c/J — 10 _1 , Tj/J = 10 -3 (/3j = Tj + T-^c^ij) an d the dephasing 
constants (7i/J,7d/J) = (0,0), (10 _1 ,0), (0, 10" 1 ) from the first to the third row. 



is an exponentially decaying function exp(— tjr\ r = tr. Another relevant temporal 
characteristic is the transfer time to the RC, defined by 

1 r°° 

r RC = - / dtr RC (l\p(t)\l)t, (61) 
6 J 
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(c) Dephasing at the donors 

□Z ■ u 

0.026 0.034 

Figure 6. (color online) Inverse relaxation time (r^J) -1 of the bi-pathway complex 
in figure [TJb) with symmetric couplings constants J = J12 = J13, parameterized by 
energy- level mismatches. In each panel, horizontal and vertical axes represent energy- 
level mismatches AQ12/J and AO13/J, respectively. The parameters are chosen as 
Frc/J = 10 _1 , Tj/J = 10 -3 (Pj = Tj + Trc^Ij) an d f ne dephasing constants 
(7i/J,7d/J) = (0,0), (HT 1 ,*)), (0, lO" 1 ) from the first to the third row. 

where (Tftc/e) (1| p(t) |1) is the probability density function of time t that an exciton 
is transferred to the RC at t. When the time dependence of the transfer to the RC is 
given by a function e(l — exp(— t/r)), r = trc- 

For a bi-pathway complex with couplings J = J12 = J13, the inverse transfer time 
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(a) Transfer Efficiency (b) Temporal Characteristics 

Figure 7. (color online) Transfer efficiency and temporal characteristics of the 
bi-pathway complex used in figure [5] and figure [6] for anti-symmetric energy-level 
mismatches AQ12/J = — AQ13/J = AQ a /J with no dephasing noise. As the 
asymmetric degree of energy-level mismatches AVt a /J is increased up to 0.1, (a) the 
transfer efficiency e increases slightly higher than 1/2 but (b) the inverse transfer 
time to the RC (trcJ) -1 decreases close to the inverse relaxation time (t#J) _1 . For 
AO a /J > 0.1, the difference between two temporal characteristics becomes negligible 
as the relaxation is mainly caused by the transfer to the RC. 




10~ 7 10" 5 0.001 0.1 10 1000 10~ 7 10~ 5 0.001 0.1 10 1000 



(a) Transfer Efficiency (b) Temporal Characteristics 

Figure 8. (color online) Transfer efficiency and temporal characteristics of the 
bi-pathway complex used in figure [5] and figure [6] for no energy-level mismatches 
A^i2 = A^i3 = with dephasing noise 7^ at the donors. No dephasing noise at 
the receiver is assumed, 71 = 0. 



to the RC (trc</) -1 and relaxation time (t#J) _1 , normalized by J, are presented in 
figure [5] and figure |6| respectively. Note that the temporal characteristics have peaks at 
the points where the transfer efficiency e is high in figure |3](a)-(c). In the first row of 
panels in figure [5] where no dephasing is assumed, the inverse transfer time is maximal 
at no energy-level mismatches and has a sharp peak along the symmetric energy-level 
mismatches AQ — A^i 2 = Afii 3 . This is due to the fast transfer to the RC of the non- 
destructive components of ^=2,3) in equations (38)-(39) (see figure The fast transfer 
occurs in a probability less than 1/2, while the destructive component of \vi) results 
in the exciton-loss at the donors in a probability of 1/2. The probability amplitudes 
of 1^=2,3) decay with rate constants \(f3\ + fid) + (— l) p=2 ' 3 Re[y] where Re[Y] is the 
real part of Y in equation (43). The inverse transfer time to the RC is determined by 
the probability of an exciton being at the receiver and for the symmetric energy-level 
mismatches it is upper bounded as 

_t & + & f 1 Q3i + /3 rf ) 4 + 4(A - /3 d ) 2 Af2 2 \ - 1 ft + fa 
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The inverse transfer time approaches the upper bound if fii, fid <C J12, J13 and AQ = 0, 
as in the case of the bi-pathway complex in figure [5J (trc</) -1 ~ + fid)/ J = 0.051. 
On the other hand, the probability amplitude of decays with a rate constant \fid 
(see equation ([42]) ) . This causes the slow relaxation of the symmetric energy-level 
mismatches (see the first row of panels in figure [6| : 



Tr 1 = 2fa 1 + 



MJ?2 + J?3)+0lfoKPl+fa) 2 +WlfaMl 2 J - 2/3d ' (63) 

For slightly-asymmetric energy-level mismatches, partially destructive interference 
occurs at the receiver molecule, which increases the transfer efficiency e slightly higher 
than 1/2 but decreases the inverse transfer time (trc</) -1 instead (see figure [7] and 
table[T]). As the asymmetric degree of energy-level mismatches increases, the interference 
effect eventually disappears and the temporal characteristics are optimized at the anti- 
symmetric energy-level mismatches A^i 2 /J = — A^i 3 /J w ±1. As in the third rows of 
panels in figure [5] and figure [6| the interference effect also disappears as increasing the 
dephasing 7^ at the donors (see figure [8] and table [T]) . 



5. Multi-Pathway Complex 

We now briefly discuss the generalization to a multi-pathway complex consisting of 
sub-complexes with multi-level antenna molecules. Suppose that a receiver is coupled 
to independent sub-complexes. When decoupled, each sub-complex yields a set of 
energy eigenvalues. It is found that if any pair of sub-complexes share a common 
energy eigenvalue, there arises destructive interference at the shared molecule, which 
funnels the excitation energies from the sub- complexes. In particular, consider the FMO 
complex consisting of seven BChl molecules, approximated as a bi-pathway complex 
where a receiver (molecule 3) is coupled to two sub-complexes, one (s = 1) consisting 
of molecules 1 and 2 and the other (s = 2) of molecules 4 to 7 (Here we use the usual 
numbering of the BChls, which was chosen by Fenna and Matthews [2Tj). We present 
the energy eigenvalues of each sub-complex in figure [9](a), based on the Hamiltonian of 
a FMO monomer of Prosthecochloris aestuarii in Ref . [22j , 



H = 
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5.1 


-61.5 


-2.5 


32.7 


280 / 



(64) 



where we shifted the zero of energy by 12230 and all numbers are given in the units 
of cm -1 . As shifting the energies of the two sub-complexes by AQ s= i j2 , we present the 
transfer efficiency without/with dephasing in figure [9](b) /(c). Here, we assume that an 
exciton is created at molecule 1. It is clearly shown by seemingly four diagonal lines in 
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(a) FMO complex (b) 7 = (c) 7 = 10 2 /188 



0.5 1 

Figure 9. (color online) (a) Energy eigenvalues of two sub-complexes of the FMO 
complex (blue for one sub-complex with molecules 1 and 2 and purple for the other 
with molecules 4 to 7) and the excited energy of the receiver (brown). Here, we shifted 
the zero of energy by 12230 (all numbers are given in the units of cm -1 ). Transfer 
efficiency of the energy-shifted FMO complex with (b) no dephasing noise and (c) 
uniform dephasing constants jj = 7, Mj. In each panel, horizontal and vertical axes 
represent the energy shifts AQ S of the sub-complexes s = 1,2, respectively, and the 
parameters are chosen as Trc = 10 3 /188 and Tj = 1/188, Vj (/3j = Tj + Trc^i?)- 

figure [9](b) that destructive interference at the receiver results in the suppression of the 
transfer efficiency. The actual suppressing lines are eight but they are overlapped due 
to the almost equal energy-eigenvalue differences of the sub- complexes. In figure [9](c), it 
is shown that dephasing noise reduces the destructive interference, which consequently 
enhances the transfer efficiency of the FMO complex (AQ S = 0) from about 77% to 96% 
when 7 = 10 2 /188. We would note that there can arise destructive interference inside a 
sub-complex if it has a multi-pathway structure depending on the excited energies and 
electronic coupling constants of the constituent molecules. 

6. Remark 

In summary, we investigated the role of energy-level mismatches between sub-complexes 
in a multi-pathway complex of a light-harvesting complex. For a single-pathway 
complex, we showed that the presence of energy-level mismatches never enhances the 
transfer efficiency, as it is likely to cause the Anderson localization that reduces the 
energy transfer through the energy pathway. However, we demonstrated that for a 
multi-pathway complex, the absence of energy-level mismatches causes the quantum 
destructive interference at the shared molecule, which funnels the excitation energies 
from the sub-complexes. We showed that the undesired destructive interference 
diminishes as introducing energy-level mismatches so that the transfer efficiency is 
significantly improved, even though they may cause the Anderson localization effect 
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and suppress the energy transfer in each energy pathway. Due to the competition 
between the localization effect and the destructive interference at the shared molecule, 
it is crucial to find an optimal amount of energy-level mismatches so as to maximize the 
transfer efficiency as minimizing both the undesired effects. We also showed that the 
quantum dephasing noise relaxes the energy localization and destroys the destructive 
interference at the shared molecule, implying that the cooperation of the energy-level 
mismatches and the dephasing noise significantly increases the transfer efficiency after 
all. This mechanism will be qualitatively unaltered even if considering non-Markovian 
noise P, [18] , as long as the role of the noise is to destroy the destructive interference 
at the shared molecule. The results presented here suggest that the energy-level 
mismatches typically given in natural photosynthetic complexes can play a crucial role 
in the efficient energy transfer. It is a timely question to investigate the optimal energy 
landscape with energy-level mismatches and non-uniform electronic couplings in a more 
elaborate dephasing model which accounts the thermal fluctuations of surrounding 
proteins and the intra/inter-molecular vibrations [23]. 
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